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Summary 

Remcom’ s XFDTD software was used to model the 
properties of a two-dimensional left-handed 
metamaterial (LHM) flat lens. The focusing capability 
and attenuation of the material were examined. The 
results showed strong agreement with both experimental 
results and theoretical predictions of focusing effects and 
focal length. Moreover, the inherent attenuation in the 
model corresponds well with the experimental results 
and implies that the code does a reasonably accurate job 
of modeling the actual metamaterial. 

Introduction 

Recently there has been considerable interest in left- 
handed metamaterials (LHM), which are also known as 
double negative metamaterials or negative index 
materials. These materials are engineered so that their 
effective electrical permittivity e and magnetic 
permeability p values are negative over a specific 
frequency range. In 1968, Veselago (ref. 1) postulated 
that a flat piece of material with simultaneous negative 
values of e and p would alter the phase and group 
velocities of electromagnetic radiation so that they 
would propagate in opposite directions and also invert 
Snell’s law of refraction, the Doppler effect, and 
Chere nk ov radiation. Pendry et al. (refs. 2 and 3) first 
showed that a three-dimensional array of intersecting 
thin wires could be used to create negative values of 8 
and that a periodic array of copper split-ring resonators 
could produce an effective negative p. In 2000 Smith et 
al. were the first to successfully combine the two arrays 
and produce a LHM which had negative values of s and 
p for a band of frequencies in the GHz range (ref. 4). 

Wilson and Schwartz experimentally tested a flat lens 
constructed of LHM (ref. 5). They created a 20- by 20- 
by 3 -unit cell array comprising copper split-ring 


resonators and thin copper wires etched onto a fiberglass 
substrate. The entire block measured 100 by 100 by 10 
mm, and a partial view is shown in figure 1 . The LHM 
block was placed in a parallel-plate waveguide as shown 
in the photograph of its experimental setup (fig. 2). On 
one side of the LHM, a coaxial cable was attached to the 
top plate of the waveguide and the center conductor of 
the cable was fed through to the bottom plate. 
Cylindrical waves with a vertically directed electric field 
radiated from this source point and passed through the 
LHM. In order to detect the microwave power 
distribution after the waves had passed through the 
LHM, a detecting probe was attached to a computer- 
controlled x-y positioning stage and was scanned over a 
distance of 100 mm in the x-direction (toward and away 
from the LHM) and over a distance of 60 mm in the y- 
direction (from edge to edge in the waveguide). A top 
view of the experimental setup is shown schematically in 
figure 3. The LHM was placed at nine locations away 
from the source so that the front focal length (designated 
as ffl in fig. 3) ranged from 30 to 70 mm. The signal was 
swept from 9.50 to 12.50 GHz in 0.10-GHz intervals at 
each location. The back focal length (bfl), which for this 
experiment is defined as the distance from the back edge 
of the block to the power maximum along the x-axis (see 
fig. 3), was measured for each frequency. The results of 
the experiment at 10.30 GHz are shown in figure 4. A 
focusing effect can be seen around the y = 0 axis. In 
addition, we find that the bfl decreases as the distance 
from the source increases. 

Modeling the geometry of the split-ring resonators and 
wires of the left-handed metamaterial flat lens would be 
too memory intensive to simulate, so the Drude model 
for explaining conductivity in materials (refs. 6 and 7) 
was used to create homogeneous frequency-dependent 
materials with negative values of both e and p. In this 
report, the effectiveness of XFDTD in simulating LHMs 
will be examined along with the properties of a variable 
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Figure 1. — Partial view of LHM block used by Wilson and 
Schwartz (ref. 5). 



Figure 2. — Experimental setup of Wilson and Schwartz (ref. 5). 
On left side of LHM block, a coaxial cable is attached to top 
plate of waveguide and the center conductor of cable is fed 
through to bottom plate. Cylindrical waves with vertically 
directed electric field radiate from this source point and pass 
through LHM. To detect microwave power distribution after 
waves had passed through LHM, a probe is attached to 
computer-controlled x-y positioning stage and scans over 
100 mm in x-direction (toward and away from LHM) and over 
60 mm (from edge to edge in waveguide). 

focal length flat lens constructed from these materials. 
The results will be compared to both earlier theoretical 
predictions (ref. 1) and Wilson and Schwartz’s 
experimental findings (ref. 5). 



Figure 3. — Top view of experimental setup of Wilson and 
Schwartz (ref. 5). Source point is designated as front focal 
point (ffp); distance from this point to LHM block is front 
focal length (ffl). Shown is a ray entering LHM block at 
angle of refraction 0 r . As it passes into LHM, the ray 
experiences negative refraction, with its angle changing 
from 0 r to 0|. As it passes out of LHM it again experiences 
negative refraction, with its angle changing from 0| to 0 r . 

This enables the radiation to focus at back focal point (bfp), 
with distance from LHM designated as back focal length 
(bfl). 

Modeling Left-Handed Metamaterial Flat 
Lens with Remcom XFDTD 

Remcom XFDTD is a full wave three-dimensional 
electromagnetic solver. The terminology “full wave” 
indicates that no approximations are used in the 
calculation of electromagnetic waves. It utilizes the 
finite difference time domain method for solving 
Maxwell’s equations. This method allows the user to 
view the electric and magnetic fields in the time domain 
and the fields’ propagation over a period of time. The 
basic unit of calculation for this method is the Yee cell 
(ref. 7). This cell is a three-dimensional rectangular box 
that places electric field values on the edges of the cell 
and magnetic field values on the sides of the cell. This 
placement of the fields allows the program to quickly 
calculate field values by using shared sides and edges 
from other cells. 

Defining the Structure Geometry 

For the first step, we need to define the geometry of 
the structure to be modeled. A screenshot of the XFDTD 
geometry screen is shown in figure 5. The largest portion 
of the screen is taken up by the geometry window. This 
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Figure 4. — Results from Wilson and Schwartz experiment 
at 10.30 GHz (ref. 5). (a) Front focal length (ffl) is 30 mm 
and corresponding back focal length (bfl) is 87.5 mm. 

(b) ffl is 50 mm and bfl is 60 mm. (c) ffl is 70 mm and bfl 
is 40 mm. 

is where the actual drawing is shown. Three different 
planes — xy, yz, and zx — are shown in small screens to 
the left as well as the three-dimensional view. The large 
geometry can be oriented to any of these views by 
simply clicking on the smaller window. The object tree 
shows each element of the geometry that has been 
created. Clicking on an element will highlight it in the 
geometry window. The shape tabs allow the user to 
create a variety of three-dimensional volumes by 
specifying properties such as width or depth as well as 
assigning predefined material properties. A two- 
dimensional editor can be used for creating complex 
planar objects and can be opened in one of three ways: 
(1) pressing the 2D Editor icon on the geometry window 
toolbar while in solid mode, (2) right-clicking in the 



Figure 5. — Screenshot of XFDTD main geometry tab, where 

model and its materials are defined. 

main view and choosing “2D Editor,” or (3) selecting 
“2D Editor” under [Tools] from the geometry window 
menu bar. This editor uses predefined shapes such as 
rectangles and ellipses, but also allows the user to draw a 
shape freehand. The drop-down menus above the buttons 
contain several features that are used to manipulate 
individual elements, although most of these features are 
also available by right-clicking on a specific part in the 
object tree. The material area is where the user is able to 
create a new type of material by defining certain 
parameters. Three default materials are defined: free 
space, perfect electric conductor (PEC), and perfect 
magnetic conductor (PMC); all other materials must be 
user defined. In XFDTD, a single material cannot be 
assigned both electric and magnetic properties. Thus to 
create a volume with both types of properties, both an 
electric and a magnetic material are assigned to occupy 
the same space in the geometry. Finally, the viewing 
tools allow the user to alter the view of the geometry by 
rotating the axes, panning, and zooming. 

It is recommended to first define the materials before 
creating the structure. This way materials can be easily 
applied to objects as they are created. In the “Material” 
area, first select either “Electric” or “Magnetic” material. 
Both options offer the user several choices to customize 
a material. The user can define a name for a material 
along with material parameters and its viewing color. 

Next the user must construct the structure of the 
geometry. This is done by using the shape tabs. XFDTD 
allows the user to create simple three-dimensional 
shapes such as rectangular boxes, spheres, pyramids and 
wires. Each object has a set of coordinates that must be 
defined before they are created. For example, a 
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rectangular box has two corner points that define its 
length, width, and depth. More complex three- 
dimensional shapes are created using Boolean 
expressions such as addition, subtraction, or intersection. 
A two-dimensional editor is available to create flat 
polygons or circles to be used in simulations. 
Unfortunately, no extrude option is available for these 
two-dimensional shapes. A computer-aided design (CAD) 
importer allows the user to bring in previously designed 
geometries to use with XFDTD and is recommended for 
more intricate geometries. XFDTD takes care of the 
meshing of the geometry once it has been imported. 


Material Characterization 

The Drude model for explaining conductivity in 
materials, first proposed by Paul Drude in 1900 
(refs. 6 and 7), was used to create the negative values of 
s and p in XFDTD. This model describes the interactions 
electrons have amongst each other and with other atoms 
(ref. 8). The materials used were created using this 
model in XFDTD, and the equations for it are 

e r (co) = e 00+ k^+-2- (1) 

1 + ycoto ycot 0 


Mesh 

We need to mesh the structure after defining its 
geometry. This involves breaking the structure into 
thousands of Yee cells. XFDTD has a default setting that 
balances memory usage with reasonable simulation time, 
which is good for most applications. It can also be set to 
mesh the structure into a user-defined number of cells, or 
the user can specify the dimensions of the cells to create 
more rectangular cell shapes. The user also has the 
option to define padding cells, or free space, surrounding 
the structure. The user can choose between a uniform 
amount of padding on each side or explicitly defining 
padding amounts on all six sides. XFDTD offers a color- 
coded view of each element and the padding surrounding 
it. Also, along the top of the view is the size of the mesh. 
It is necessary to consider this because the size will be a 
large factor in determining the run time of the 
simulation. The user must be aware of the amount of 
RAM the computer has because meshes that are too 
large will not run at all. 

XFDTD comes with an adaptive mesh feature. Here, 
the user is able to specify certain regions of the mesh to 
refine. It is possible to define either exact cell 
dimensions or the ratio to the main mesh cell size in the 
specified direction. For example, if the main cell size is 
1 mm, a specified ratio of 4:1 in the z-direction will 
provide a z-dimension length of 0.25 mm for the cells in 
the refined region. XFDTD automatically adjusts the 
cells sizes at the boundaries of the adaptive regions to 
transition back to the main cell size. 

After the mesh is created the user can alter the mesh 
by using the [Mesh] edit tool. This allows the user to 
assign a different material to specific cells. The mesh 
edit feature is useful if a small change is needed, but 
large changes should be made by creating a new object. 




1 + y'roto 7®h) 


( 2 ) 


with Soo and p* being the infinite frequency relative 
permittivity and permeability, respectively; s s and p s 
being the static relative permittivity and permeability, 
respectively; a and c,„ being the electrical and magnetic 
conductivity, respectively; co being the angular 
frequency; and t 0 being the relaxation time. To calculate 
the refractive index, first rewrite the equations into their 
real and imaginary parts: 


8 r ({o) = E , -y8" 
Pr(®) = P'-7F" 


0) 

(4) 


such that 


s 


= s„ 


(£, -O 

1 + (ofr 0 )- 


„ _ ( £ s ~ £ °o) + g 
1 + (Otfro ) CDSq 


(5) 

( 6 ) 


The same form can be applied to the magnetic case. The 
refractive index n is defined by equation (7): 


n( co) = ± i ys(co)p(co) 


(7) 


Electromagnetic waves will not propagate through a 
material whose refractive index is imaginary. Therefore, 
the signs of e(co) and p(co) must be the same. A negative 
sign for the square root must be chosen when both s and 
p are negative (ref. 9). By substituting equations (3) and 
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(4) into equation (7), the refractive index is 

n = ±V(s'-7V')0i'-yy') (8) 

After solving for the real part of n in equation (8), the 
following equation is obtained 


2 

*^(sy)-(s"p") + V(s'Vs" 2 )(p ,2 +p" 2 7 (9) 

where n' is the real part of the index of refraction. It can 
be seen that with proper selection of values, the real 
portions of the permittivity and permeability can become 
negative. These numbers cannot be selected randomly in 
order to achieve negative values because of software 
calculation limitations. These limitations are that the 
time step of the simulation must be smaller than 
3 percent of the relaxation time, e, and g-, should not be 
much larger than 1 , and finally, if e v < b v and p. s < |i x — 
which is necessary to achieve negative s(co) and p(co) — 
then 



x So S oo ) 
a > 

U 


and 


a 


m 


> Po(lU - Too) 
f o 


( 10 ) 


( 11 ) 


To calculate acceptable values, this information was 
entered into a spreadsheet, and the values were altered 
until all conditions were met. The resulting values were 
Boo = 1.001, geo = 1.001, b s = -4269, |i s = -4269, a = 
54.02, a m = 7.667xl0 6 , and to = 7.0xlCT 10 . Using these 
values, a graph was created to show how the values of n' 
vary over a very narrow band of frequencies. This is 
shown in figure 6. 


Signals and Ports 

A port is any side of a Yee cell in the mesh where a 
signal will be applied or a measurement will be taken. 
These ports can be added in two ways. The first is by 
simply going under the [Run Parameters — » Component 
Ports] tab and entering in the cell location, orientation, 
polarity, and component type. A second and easier way 



Figure 6. — Real part of index of refraction n' as function of 
frequency calculated with Drude model. 


is to view the geometry in the [Mesh] mode and right- 
click on the desired cell. Next select “Edit Port” on the 
menu that appears. This takes the user to the 
[Component Ports] tab as well, but all of the cell 
information has already been entered. 

After all port parameters have been set, the type of 
signal must be defined. This is done under the [Run 
Parameters — » Source Type] tab. Here, the user can 
define a plane wave, discrete source, or a Gaussian beam 
and the polarity of the signal. Next, the shape of the 
signal is selected under the [Waveform] tab. The user 
can choose to examine a broadband result with a 
Gaussian pulse, Gaussian derivative, or modulated 
Gaussian or look at a steady-state result with a sinusoid. 
The user also has the option to create a unique type of 
signal by using the “User Defined” option. Finally, 
boundary conditions can be set under the [Outer 
Boundary] tab. A boundary can be a Perfect Electric 
Conductor (PEC), Perfect Magnetic Conductor (PMC), 
Liao, or Perfectly Matched Layer (PML). The PEC 
boundary is used when there is a conducting surface at 
the boundary. Both PEC and PMC boundaries can be 
used to take advantage of geometry symmetries to 
reduce the size of the mesh and the calculation time (ref. 
10). The Liao boundary (ref. 11) is an absorbing 
boundary, which is usually quick and accurate, but it 
does produce instabilities in the calculations if objects 
are placed too close to the boundary. The PML boundary 
(ref. 12) is an absorbing boundary that is more robust 
than the Liao boundary, but the method of calculation is 
more complex and adds a significant amount of time to 
the simulation. 

Request Results 

In the final step before the simulation can be run, the 
user must decide in which plane(s) to view the output 
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data. Once again, this can be done in two ways. Under 
the [Request Results] tab, the user selects the type of 
output to be calculated. Common types are near-zone 
fields and far-zone fields, but the transient field 
snapshots option is also very useful. This option saves a 
sequence of the electric field E, magnetic field H, 
current density J, and Poynting vector S values that can 
be viewed individually or as an animation in the 
[Geometry — » View] tab after calculations. 

Postcalculation 

After calculations have been run, there are a variety of 
results that can be examined. Already mentioned are the 
movies of the fields, but they are made even more useful 
by providing the feature to view them in vector format. 
This allows the user to view not only the magnitude of 
the field but also its orientation as it rotates through 
space. Also useful are the plots that XFDTD offers. The 
user has several options to view such as voltage versus 
time or impedance versus frequency. There are also 
options for far field calculations or fast Fourier 
transforms. S-parameter calculations are available for 
examining multiple port systems, but these must be 
turned on before calculations are run. 





(a) 




Results 

A block of electric material and a block of magnetic 
material were created so that both blocks occupied the 
same space. A sinusoidal source was added on the left 
side of the block with a frequency of 10.5 GFlz. Only the 
transient fields were requested for these simulations. 
Figure 7 shows the results from simulations with the 
metamaterial blocks 30, 50, and 70 mm to the right of 
the source. The material is represented by the white 
outline, and the wave is traveling left to right. When the 
waves travel across the boundary of the material, the 
shape of the resultant waves appear to be the mirror 
image of the incident wave, and a focusing effect can be 
seen inside the material. The sharpness of the focal 
point, though, has been slightly degraded. In the 
animation of these waves, they actually appear to be 
traveling in opposite directions; the reason for this is the 
opposite directions of the group and phase velocities 
created by the negative values of s and p (ref. 1). After 
the wave comes to a focus, it once again reverts back to 
the shape of the original wave, although the strength of 
the wave has attenuated. When the wave encounters the 
second boundary, the shape is once again reversed and 
begins to focus. This focal point also produces a highly 


Figure 7. — Results of XDTFD simulations for LFIM flat lens, 
with LHM block to right of source. A focal point, inside 
LFIM block, is indicated by ■. Back focal point, on right 
side of LHM block, is indicated by A. Intensity of field is 
represented by color, with red being the strongest and 
blue the weakest, (a) LHM block is 30 mm from source 
and back focal length (bfl) is 58 mm. (b) LHM block is 
50 mm from source and bfl is 40 mm. (c) LHM block is 
70 mm from source and bfl is 27 mm. 

attenuated wave, and it can clearly be seen that the 
resolution of the focal point has greatly deteriorated. 

A flat lens made from a LHM shoidd theoretically be 
able to change its focal length by adjusting the distance 
of the material from the source (ref. 5). This was tested 
by running three simulations with the LHM lens 30, 50, 
and 70 mm away from the source. The results for the 
three tests are shown in figure 7. These results 
correspond to the theoretical and experimental results 
obtained by Wilson and Schwartz (ref. 5) that were 
shown in figure 4. It can be seen that as the lens moves 
away from the source, the back focal length (bfl) 
becomes shorter. With the lens 30, 50, and 70 mm from 
the source, the bfls were approximated to be 58, 40, and 
27 mm, respectively. It can also be noticed that the 
second focal point becomes less defined as the LHM 
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Figure 8. — LHM flat lens focal lengths modeled with XTDFD 
compared with theoretical and experimental results of 
Wilson and Schwartz (ref. 5). 


block is moved further away from the source. Figure 8 
shows these results along with Wilson and Schwartz’s 
results. Both the XFDTD and experimental results agree 
moderately well with the theoretical values. 

A large amount of attenuation, as much as 16.1 dB, 
was also noticed in the experimental results. A higher 
degree of attenuation of about 25 dB was noticed in the 
XFDTD simulations. The extent of the attenuation of the 
signal in XFDTD is determined by the parameters in the 
Drude equations. Less attenuation and better matching 
with experiment might be possible by using different 
values in the Drude equations with frequencies outside 
of the range considered. 

Conclusions 

A variable focal length flat lens constructed from 
LFIMs was simulated in XFDTD (Remcom). Drude 
equations can be used to simulate negative values of 
electrical permittivity and permeability if the parameters 
are set very carefully. XFDTD exhibited the ability to 
calculate the fields for LHMs and illustrated the effect of 
a wave having opposite values of phase velocity and 
group velocity. A flat lens made from this material did 
indeed show a focusing effect despite very high 
attenuation in all simulations. The focal length was 
varied by adjusting the distance of the lens from the 
source. At 30 mm from the source the back focal length 


(bfl) was approximately 58 mm, at 50 mm from the 
source bfl was 40 mm, and at 70 mm from the source the 
bfl was 27 mm; these compared moderately well with 
theoretical values. 
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